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^ Abstract: An analysis of transition from chaotic to nonchaotic behavior and synchronization in 
^ an ensemble of systems driven by identical random forces is presented. The synchronization phe- 
^ homenon is investigated in the ensemble of particles moving with friction in the time-dependent 
Q\ potential and driven by the identical noise. The threshold values of the parameters for transition 
■^ from chaotic to nonchaotic behavior are obtained and dependencies of the Lyapunov exponents and 

,.^^ower spectral density of the current of the ensemble of particles on the nonlinearity of the systems 
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O and intensity of the driven force are analyzed, 
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H , A shortened version of the paper: 



SYNCHRONIZATION BY IDENTICAL NOISE 



1. Introduction 

Often trajectories of tlie nonlinear dynamical systems are very sensitive to the initial conditions and 
unpredictable, i.e. the systems are chaotic. These systems exhibit an apparent random behavior. It 
might be expected that when turning on additional random forces make their behavior even "more 
chaotic" . However, a transition from chaotic to nonchaotic behavior in an ensemble of particles with 
different initial conditions bounded in a fixed external potential and driven by an identical sequence 
of random forces was observed by Fahy & Hamann [1992] and recently analyzed theoretically and 
numerically by Kaulakys & Vektaris [1995a,b] and Chen [1996]. It has been shown that the ensemble 
of trajectories in such a case may become identical at long times. The system becomes not chaotic: 
the trajectories are independent on the initial conditions. The similar effects have been observed in 
the different systems as well [Yu et ai, 1990 and Maritan & Banavar, 1994a] and have resulted to 
same discussion concerning the origin and causality of such nonchaotic behavior [Pikovsky, 1994, 
Maritan & Banavar, 1994b and Gade & Basu, 1996]. 

Moreover, the observed effect resembles a phase transition but does not depend crucially on the 
dimension of the space in which the particles move. This phenomenon has some importance for 
Monte Carlo simulations and can influence on the clustering of particles process. 

It should be noted that Maritan & Banavar [1994a] have analyzed the similar effect using the 
Langevin equation, however, in the limit in which the time separation between successive forces 
become small compared to any characteristic macroscopic time of the system and two logistic maps 
linked with a common noise term. 

Here we analyze the similar phenomenon in the ensemble of particles moving with friction in the 
time-dependent potential and driven by the discrete identical noise. We define the threshold values 
of the parameters for transition from chaotic to nonchaotic behavior and investigate dependencies 
of the Lyapunov exponents and power spectral density on the nonlinearity of the systems and 
character of the driven force. Our theoretical analysis is based on the mapping form of equations 
of motion for the distance between the particles and the difference of the velocity of the particles 
while numerical calculations are performed according to the derived mapping equations as well as 
directly calculating the system's trajectories and the Lyapunov exponents. The mapping analysis 
results in the conclusions very close to those obtained from the direct simulations and numerical 
calculations 



2. Models and Theory 

Consider a system of particles of mass tti moving with friction according to Newton's equations 

dr dw 1 dV (r, t) 

dt ^ dt m dr 

in the time dependent potential V{r, t), e.g. in the potential V {x, t) = x'^ — x^ — ax sin ut, and with 
the friction coefficient 7. 

At time intervals r the particles are partially stopped and their velocities are reset to the mixture 
of some part a of the old velocities with the random velocity v™": v"^"" = av°'^ + v™", where k 
is the stop number. Note that v™" depends on the stop number k but not on the particle. The 
simplest and most natural way is to choose the random values of velocity v™" from a Maxwell 
distribution with ksT = m = 1, i.e. from the Gaussian distribution of variance a^ = 1. Figure 1 
illustrates the difference of evolution of the ensemble of particles with randomly distributed (from 
the Gaussian distribution of variance o"^ = 1) initial conditions and perturbed by the replacement 
^new _ (y^^oid _|_ ^ran j-^j, different values of the time interval r between such perturbations. We see 
the transition to one (common for all particles) trajectory for sufficiently small time interval r. 

Theoretically a transition from chaotic to nonchaotic behavior in such a system may be detected 
from analysis of the neighboring trajectories of two particles initially at points ro and Tq with starting 
velocities Vq and Vq. The convergence of the two trajectories to the single final trajectory depends 
on the evolution with a time of the small variances Arj^ = r^ — r^ and Av^ = v^ — v^ . From formal 
solutions r = r(rfc,Vfe,t) and v = v(rfc,Vfe,t) of the Newton's equations with initial conditions 
Y = Yk and V = Vfc at t = it follows the mapping form of the equations of motion for Ar and Av 
[Kaulakys & Vektaris 1995a,b]: 

Arfc+i \ rr./ X / Arfc \ , . 

Av,,J=T("'^'='^^'^'=HavJ ^2) 



where the T matrix is of the form 



vr Lt^ vv 



/ (9r Oy \ 

dYk dvk 

(9v (9v 

V dYk dvk ' 



(3) 



and the time interval r^ may depend on the step k. 

For one-, two- and three-dimensional systems the dimension of the T matrix is 2, 4 and 6, 
respectively. 



Matrix elements Trr and T^v satisfy the equation 

d^Tr _ If d\dV (r, t) 

J- r 



dt'^ m y ^ drj dr 

and initial conditions at t = 



— (4) 

r=r(rfc,Vfc,t) 



Trri^'k, Vfc, 0) — Tvv — 1, Trv — T^r — 
Trri^k, Vfc, 0) = Tvv = 0, Trv = 1, T^r '- 



) 



(5) 



m dr^ 

while Tvr = Trr and Tw = Trv Here and further the points over the letters express the derivatives 
with respect to the time. 

Further analysis is based on the general theory of dynamics of classical systems represented as 
maps: we can calculate the eigenvalues fik of the T matrix for each step and evaluate the averaged 
Lyapunov exponents or KS entropy of the system 

/I \ 1^1 

o-fe = ( — ln|;Ufe|) = lim — V — In |/ifc(rfc, v^, r^)! (6) 

A criterion for transition to chaotic behavior is 

'-'max U . yi ) 

Comparisons of the threshold values Tc for transition to chaos according to Eqs. (2)- (7) with 
those from the direct numerical simulations indicate to the fitness and usefulness of the method (2)- 
(7) for investigation of transition from chaotic to nonchaotic behavior in randomly driven ensemble of 
systems bounded in the fixed external potential without the friction [Kaulakys & Vektaris 1995a,b]. 

3. Results of calculations 

Here we calculate the Lyapunov exponents directly from the equations of motion and linearized 
equations for the variances of coordinate and velocity. Further we extend the same analysis for the 
systems with friction in the regularly time depending external field and perturbed by the identical 
for all particles random force. 

Figures 2-5 represent an extensive analysis of the autonomous system based on the numerical 
solutions of the differential equations of motion. Figure 2 shows quite similar behavior of the 
Lyapunov exponent to that calculated from the mapping equations of motion (curve (a) of Fig. 
2 in the paper [Kaulakys & Vektaris 1995a]). Figures 3-5 represent dependences of the Lyapunov 



exponents on the different parameters of the model. Areas of the parameters for which the Lyapunov 
exponents are negative corresponds to the nonchaotic Brownian-type motion. 

In general, motion in the nonautonomous Duffing potential with friction describe equations 

V = 2x — 4x^ — 7f + a sin tot, x = v. (8) 

For a = and 7 = Eqs. (8) represent motion in the fixed external potential. As it was 
mentioned above, in the paper by Kaulakys & Vektaris [1995a] some theoretical and numerical 
analysis of this model was fulfilled on the bases of the mapping equations (2)-(7). 

Figure 6 illustrates evolution of the ensemble of particles with randomly distributed initial 
conditions in the nonautonomous Duffing potential with friction and perturbed by the replacement 
f "'^"'(/cr) = av°^'^{kT) + f ™", k = 1,2, .... We observe a transition from the actual chaotic dynamics 
for large r to the nonchaotic common for all particles trajectory with the decrease of r. In Fig. 
7 we show the dependence on r of the Lyapunov exponents for the motion in the nonautonomous 
Duffing potential with friction. For the values of parameters corresponding to the positive Lyapunov 
exponents, i.e. without the random perturbation (r -^ 00), the system is chaotic. The negative 
Lyapunov exponents for small r indicate to the nonchaotic Brownian-type motion. 

4. Spectrum of the current noise 

As it has already been observed in the paper by Kaulakys & Vektaris [1995b] such systems exhibit 
the intermittency route to chaos which provides sufficiently universal mechanism for l//-type noise 
in the nonlinear systems. Here we analyze numerically the power spectral density of the current 
of the ensemble of particles moving in the closed contour and perturbed by the common for all 
particles noise. The simplest equations of motion for such model are of the form 

V = F — 'jv, X = V (9) 

with the perturbation given by the resets of velocity of all particles after every time interval r 
according to the identical for all particles replacement v'^'^^ikr) = av°^'^{kT) + t;™", k = 1,2,.... 
For sufficiently small r we observe the current power spectral density S (/) dependence on the 
frequency / close to the 1// -dependence (Fig. 8). It should be noted that such spectral density 
dependence is nonsensitive to some additional (nonlinear) terms in the equation for velocity. The 
essential condition for the l//-type dependence of the current power spectral density is the random 
sufficiently strong perturbation of the particles' velocities (see also Kaulakys & Meskauskas [1997] 
for analysis of other systems and different perturbations). 



5. Conclusions 

From the fulfilled analysis we may conclude that, first, synchronization and transition from chaotic 
to nonchaotic behavior in ensembles of the identically perturbed by the random force nonlinear 
systems may be analyzed as from the mapping form of equations of motion for the distance between 
the particles and the difference of the velocity as well as from the direct calculations of the Lyapunov 
exponents and, second, the model of Fahy & Hamann [1992] may be generalized for the ensemble 
of particles moving with friction in the time- dependent potential. 

The transition from chaotic to nonchaotic behavior in the ensemble of particles moving with 
friction may also be observed in the more realistic case of motion without the momently stops of 
the particles. On the other hand, the motion in the time-dependent potential without the random 
periodic perturbations even in the one-dimensional case may be chaotic or nonchaotic, depending 
on the system's parameters. Therefore, such generalization of the model allows to investigate the 
synchronization phenomenon and transition from chaotic to nonchaotic behavior effect in an ensem- 
ble of systems driven by identical random forces on the base of relatively simple (one-dimensional) 
models for larger variety of the system's dynamics. 

Moreover, an ensemble of systems linked with a common external noise may exhibit the 1//- 
type fluctuations. Our model may easily be generalized for systems driven by any random forces 
or fluctuations. On the other hand, the phenomenon when an ensemble of systems is linked with a 
common external noise or fluctuating external fields is quite usual. Thus, an ensemble of systems 
in the external random field may provide a sufficiently universal mechanism of l//-noise. 
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Figure captions for the paper 

Fig. 1. Illustration of the difference of evolution in the {x,v) -space of the ensemble of particles 
in the autonomous Duffing potential V{x) = x'^ — x^ with randomly distributed initial conditions 
and perturbed at time intervals r by the replacement f"'^'"(A;r) = Q;f°'^(A;r) +f™"', A; = 1,2, ...with 
a = 0.5. For the relatively large r = 2 there is no transition to the common trajectory, for smaller 
r = 0.8 the clustering process of particles with different initial conditions is relatively slow while 
for sufficiently small r = 0.6 a collapse to the common trajectory at the time moment t = 100 is 
evident. 

Fig. 2. Lyapunov exponent (multiplied by r) from the direct calculations vs the time r between 
the resets of the velocity v"''^'^{kr) = av°^'^{kT) + /Sf^"", k = 1,2, ... for motion in the autonomous 
Duffing potential. 

Fig. 3. Lyapunov exponent A for motion in the autonomous Duffing potential vs the time interval 
r between the resets of the velocity f "^""(/cr) = av°^'^{kT) + j3v'^'^, k = 1,2, ... for different values of 
the parameter a. 

Fig. 4. As in Fig. 3 but for different values of j3. 
Fig. 5. As in Fig. 3 but for small values of j3. 

Fig. 6. As in the Fig. 1 but for motion according to Eq. (8) in the nonautonomous Duffing potential 
with 7 = 0.07, a = 5, q; = 0.8 and (3 = 1. K transition from the actual chaotic (at r = oo) to the 
nonchaotic dynamics with the decrease of the time interval r between the resets of the velocity is 
observable. 

Fig. 7. Lyapunov exponent A vs the time r between the resets of the velocity f"^'"(A;r) = av°^'^{kT) + 
/9f ™"', k = 1, 2, ... for different values of the parameter a for motion in the driven Duffing potential 
with friction according to Eq.(8). 

Fig. 8. The power spectral density of the current of the ensemble of particles moving according 
to Eq. (9) with F = 1, 7 = 0.1 and perturbed by the common for all particles noise f"'^"'(A;r) = 
av°^'^{kT) + t"™", k = 1,2, ... with a = 1 and different values of r. The dense lines represent the 
averaged spectra. 
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